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VELOCITY  OP  DETONATION  PROM  STREAK  CAMERA  RECORDS 

By 

John  0.  Erkman 

ABSTRACT:  Some  of  the  problems  associated  with  measuring  the 
detonation  velocities  of  cylindrical  charges  are  discussed.  The 
streak  camera  is  used  as  the  recording  instrument  so  that  data  re¬ 
duction  requires  differentiation  of  numerical  data.  The  emphasis  Is 
on  methods  of  determining  if  the  records  represent  steady  or  nonsteady 
events. 


Carl  Boyars,  Chief 
Advanced  Chemistry  Division 
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U.  S.  NAVAL  ORDNANCE  LABORATORY 
White  Oak,  Silver  Spring,  Maryland 


i 


NOLTR  68-117  19  September  1968 

VELOCITY  OP  DETONATION  PROM  STREAK  CAMERA  RECORDS 


This  work  was  carried  out  under  ORDTASK  033  102  F009  06  01  which  is 
a  systematic  investigation  of  the  explosive  behavior  of  composite 
propellant  models.  The  report  will  interest  those  who  have  to 
measure  and  correct  detonation  velocities  of  exploding  charges  by 
use  of  streak  cameras. 
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VELOCITY  OP  DETONATION  PROM  STREAK  CAMERA  RECORDS 

INTRODUCTION 

I.  The  need  for  accurate  values  of  the  detonation  velocity  in  the 
study  of  organic  high  explosives  is  discussed  by  Priced  The  detona¬ 
tion  velocity  is  an  important  parameter  in  the  study  of  inorganic 
materials  such  as  ammonium  perchlorate  for  the  same  reasons.  Methods 

for  measuring  the  velocity  of  detonation  are  reviewed  by  Campbell  et 
o 

al  with  special  emphasis  on  chronographic  methods.  In  a  continuing 
program  at  the  Naval  Ordnance  Laboratory,  the  detonation  velocities 
of  both  organic  and  inorganic  materials  is  determined  by  the  use  of  a 
streak  camera.  This  report  is  primarily  concerned  with  the  problem  of 
extracting  as  much  information  as  possible  from  the  streak  camera  re¬ 
cord. 

The  principal  material  being  studied  in  the  present  project  is 
ammonium  perchlorate  (AP)  loaded  in  cylindrical  charges  at  densities 
ranging  from  0.5  to  1.55  g/cc.  Mixtures  with  AP  are  also  studied, 
the  second  component  being  wax,  aluminum  powder,  or  an  organic  high 
explosive.  Other  materials  included  in  the  program  are  nitroguanidine 
and  dinitrotoluene.  The  study  of  such  a  variety  of  materials  at  vary¬ 
ing  densities  presents  many  challenging  problems  In  the  preparation  of 
charges  and  In  the  use  of  the  streak  camera,  however,  preparation  of 
the  charges  and  the  handling  of  the  streak  camera  will  not  be 
discussed. 

The  most  challenging  aspect  of  the  reduction  of  the  data  from 
the  streak  camera  in  the  program  Is  that  of  determining  If  the  detona¬ 
tion  was  steady.  It  is  relatively  simple  to  extract  a  velocity  from 
the  data.  Steadiness  must  be  determined  from  the  acceleration  of  the 
event.  This  requires,  in  effect,  that  the  numerical  data  be  dif¬ 
ferentiated  twice.  Various  methods  are  discussed  which  are  being 
used  in  classifying  the  records  as  those  of  steady  or  nonsteady 
processes . 

Also  discussed  in  this  report  are  the  effects  of  overboostering 
the  acceptor  charge  and  the  effect  of  the  geometry  of  the  charge  on 
the  results.  The  latter  comes  about  because  the  streak  camera  records 
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the  progress  of  the  detonation  along  the  exterior  of  the  cylindrical 
charges .  Limitations  of  the  firing  facility  restrict  the  length  to 
diameter  ratio  (L/d)  to  about  3 <>5  for  3-0  inch  diameter  charges.  For 
such  charges,  the  velocities  obtained  by  observing  the  periphery  of 
the  charges  must  be  corrected  to  give  the  actual  velocity  which  would 
be  observed  if  the  axis  of  tne  charge  could  be  viewed t 

II.  CURVATURE  OF  A  STREAK  GAMER/.  RECORD.  As  noted  above,  the  nature 
of  the  research  which  prompted  this  discussion  requires  the  accurate 
determination  of  the  velocities  of  steady  detonations  in  cylindrical 
charges  of  propellants.  Unfortunately,  there  is  no  a  priori  knowledge 
concerning  the  steadiness  of  the  reaction.  This  knowledge  must  come 
from  an  examination  of  the  streak  camera  re  cord  which  can  be  done  by 
comparing  the  trace  with  a  straight  edge.  But  how  much  can  the  trace 
deviate  from  the  straight  edge?  As  a  start  in  answering  the  question, 
it  is  helpful  tc  calculate  the  radius  of  curvature  of  a  typical  trace. 
The  radius  of  curvature,  R,  is  given  by 

R  »  [1  +  U' ) 2  3 3/2  (1) 

--  x1 

where  x‘  and  x"  are  the  first  and  second  derivatives  of  a  function  of 
time,  t,  which  is  fitted  to  the  record.  For  simplicity,  assume  that 
the  trace  can  be  described  by 

2 

x  =  a  +  bt  +  ct  , 

so  that  the  derivatives  are  easily  calculated.  The  situation  is  dia¬ 
grammed  in  Fig.  1  where  the  trace  is  represented  by  the  curve.  The 
straight  line  (a  chord)  passes  through  the  beginning  and  end  of  the 
trace.  If  velocity  is  vQ  at  (0,0)  and  v-^  at  (t][,  x^),  then 

vr.  -  v^  b  -  (b  +  2c1^)  ~2ct^ 

=  P  =  b  =  ~b 

or 

c  =  -bP  /  Pt1,  (3) 

so  that  c  can  be  evaluated  for  typical  values  of  t.^  and  b  for  ap¬ 
propriate  values  of  P.  Note  that  b  is  the  velocity  at  (0,0)  in  Fig.  1. 
Now  the  equation  for  the  radius  cf  curvature  becomes. 
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R  ”  |l  +  b2  (1  -  P/2)2 1  3,,2(  (it) 

where  the  first  derivative  was  evaluated  at  t^/2.  When  the  change 
In  velocity  Is  zero  (P  »  0),  the  radius  of  curvature  is  infinite. 

Curves  showing  the  relation  between  R  and  P  (in  percent)  are  given  in 
Pig.  2.  For  two  cases,  the  trace  is  assumed  to  be  inclined  at  an 
angle  of  45°  to  the  time  axis.  For  one  of  these  traces,  t^  is  repre¬ 
sented  by  a  length  of  5.0  cm,  for  the  other,  t1  is  2.5  cm.  The  third 
curve  is  for  0  »  60°,  and  t^  -  2.9  cm.  These  values  are  fairly  typi¬ 
cal,  that  is,  the  trace  for  many  experiments  is  inclined  at  45°  and  is 
5.0  cm  long  in  the  time  direction.  Note  that  t^  is  given  in  units  of 
length  so  that  the  results  can  be  applied  directly  to  the  film  record. 
In  these  units,  the  value  of  b  is  simply  tan  0  with  no  conversion 
factor.  The  value  of  R  for  P  =  -2#  ranges  from  548  to  696  cm. 

Values  of  R  are  not  as  helpful  as  are  values  of  h,  see  Fig.  1. 

The  distance  from  the  circle  to  the  straight  line  is 

h  *  R  -  0.5  ^ 4R2  -  L2  (5) 

where  L  is  the  length  of  the  chord  of  the  trace.  Figure  5  shows  how 
h  varies  with  P  for  the  three  cases  described  above.  For  P  «  2#,  the 
values  of  h  is  about  -0.011  cm  for  0  =  45°,  t^  =  5  cm.  This  means 
that  in  order  to  detect  a  2#  decrease  in  velocity,  one  must  detect  a 
separation  of  0.01  cm  between  the  straight  edge  and  the  trace.  For 
ideal  traces  this  may  be  possible.  But  where  the  edge  of  the  trace 
is  fuzzy  or  blurred,  the  situation  is  difficult.  Magnification  in¬ 
creases  the  separation  but  may  be  of  no  real  help  for  indistinct 
records.  The  situation  is  worse  for  the  "half  record”,  i.e.,  for 
ti  =  2.5cm,  and  0  =  45°.  A  separation  of  0.004  cm  must  be  detected 
for  P  =  2#. 

Curves  such  as  those  discussed  above  help  in  the  classification 
of  records  as  representing  constant  or  varying  velocity.  After  a  pre¬ 
liminary  screening  by  direct  observations,  those  records  which  seem 
to  represent  steady  phenomena  are  investigated  further.  One  step  in 
the  analysis  is  a  determination  of  the  constants  for  the  best  quadratic 
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fit.  Once  the  constants  b  and  c  are  evaluated,  the  change  in  velocity, 
P  can  be  calculated  for  the  record.  The  result  may  be  useful  in  con¬ 
firming  the  original  classification  of  the  record. 

III.  QUALITY  OF  THE  STREAK  CAMERA  RECORD.  If  all  records  were  near 
perfect,  data  reduction  would  present  fewer  problems.  Most  records 
have  some  defects,  due  primarily  to  the  quality  of  the  charge.  For 
example,  if  the  surface  of  the  charge  is  pitted,  or  if  one  component 
in  the  mix  tends  to  clump,  the  record  will  not  be  smooth.  Also,  if 
the  charge  is  assembled  from  a  set  of  pellets,  the  record  may  not  be 
smooth  because  each  pellet  detonates  ar.  a  .slightly  different  velocity. 
The  difference  in  velocity  could  be  caused  by  variation  of  density 
from  pellet  to  pellet.  Hence  data  reduction  techniques  must  be  suf¬ 
ficiently  versatile  to  handle  records  of  different  quality. 

Data  reduction  follows  the  common  route  of  digitalizing  the  re¬ 
cord  and  performing  numerical  operations.  The  records  are  read  on  a 
projecting  machine  having  two  cross  wires.  The  inter- section  of  the 
wires  is  placed  on  the  trace,  and  a  card  is  punched.  Because  of  non¬ 
uniformities  in  many  of  the  records,  the  data  are  not  read  at  equal 
intervals  in  t.  This  is  somewhat  unfortunate,  because  many  techniques 
for  data  reduction  are  simpler  for  data  spaced  equally  in  one  of  the 
variables.  Fortunately,  more  generalized  treatments  can  be  obtained 
at  low  cost  with  modern  computing  equipment. 

One  way  to  examine  the  digital  data  for  quality  is  to  take  the 
ratio  of  the  change  of  x  to  the  change  of  t  between  each  succeeding 
pair  of  points.  This  is  a  crude  form  of  differentiation  and  is  cer¬ 
tainly  permissible  when  it  is  recalled  that  the  record  is  nearly 
straight.  The  occasional  wide  variation  of  Ax/A t  indicates  that  an 
error  has  been  made  in  reading  a  point,  or  that  the  record  is  rough. 

It  is  possible  to  set  some  criteria  on  which  bad  points  can  be  dis¬ 
carded.  This  has  not  yet  been  developed  for  the  current  program. 

After  the  data  are  on  punched  cards,  a  computer  is  used  to  de¬ 
termine  the  best  straight  line  which  fits  the  data.  Tnis  provides 
another  means  of  studying  the  scatter  of  individual  points.  If  one, 
or  a  few  points  lie  off  the  line  by  an  inordinate  amount,  the  points 
should  be  either  discarded,  or  given  a  lower  weight  in  subsequent  data 
reduction. 
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The  slope  of  the  straight  line  which  fits  the  data  Is  propor¬ 
tional  to  the  velocity  of  the  phenomenon  being  studied.  The  computer 
program  provides  an  estimate  of  the  standard  deviation  of  the  co¬ 
efficients  of  the  function  used  in  the  fit.  Hence  the  standard  de¬ 
viation  of  the  velocity  is  available.  Values  of  the  standard  devia¬ 
tion  in  the  velocity  are  being  studied  for  correlation  with  the 
occurrence  of  bad  points,  fuzzy  traces,  etc. 

IV.  CORRECTION  FOR  CHARGE  GEOMETRY  AND  BOOSTER.  The  streak  camera 
records  the  progress  of  the  detonation  along  the  exterior  of  the 
charge.  When  the  charge  length  to  diameter  ratio  is  5  or  less,  the 
velocity  derived  from  the  record  must  be  corrected  to  give  the  velocity 
along  the  axis  of  the  charge.  This  correction  has  been  discussed 
previously  by  Clairmont  and  Jaffe^.  The  situation  is  represented  by 
the  diagram  in  Fig.  4,  which  shows  the  cross-section  of  a  charge. 

Point  initiation  is  assumed  to  have  occurred  at  point  A,  following 
which  the  detonation  front  remains  spherical.  At  a  later  time  the 
front  is  represented  by  the  arc  of  the  circle,  BCD.  Using  the  point  A 
as  the  origin,  the  distance  along  the  axis  of  the  charge  to  the  in- 
stanteous  location  of  the  front  is  y.  Similarly,  using  E  as  an 
origin,  the  distance  along  the  edge  of  the  charge  to  the  front  Is 
represented  by  z.  Note  that  AE  is  the  radius  of  the  charge,  d/2,  and 
that  AC  =  AC  and 

y2  =  (AE)2  +  (DE)2  =  (d/2)2  +  z2  (6) 

Differentiating  and  substituting  for  y  gives 

y'  =  z'  /V  1.0  +  0.25  d2/y2  (7) 

where  z’  Is  the  velocity  determined  from  the  streak  camera  record. 

Let  D  =  y',  the  velocity  along  the  axis  of  the  charge  and  V  =  z',  the 
observed  velocity.  There  results 

D  =  V/  ^1  +  0.25  (d/L)2  (8) 

where  L  is  the  length  of  the  charge  up  to  the  point  observed. 

The  assumption  concerning  the  curvature  of  the  detonation  front 
has  been  discussed  previously-^  and  appears  to  be  justified  for  charges 
of  tetryl.  It  has  not  been  checked  for  the  propellants  now  being 
investigated. 
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2  -l/p 

Values  of  the  correction  factor,  F  =  (1.0  +  0.25  (d/L)  ) 
are  given  in  Table  I.  In  practice  the  combined  length  of  the  accep 
tor  and  donor  charges  is  10  inches.  For  a  5  inch  diameter  charge, 
L/d  =  5.5  and  P  =  0.989  so  the  correction  amounts  to  about  1$. 

TABLE  I 

Correction  Factor  P  as  a  Function  of  L/d 


1  0.8940 

2  0.9700 

5  0.9864 

4  0.9925 

6  0.9965 

8  0.9980 

10  0.99875 

11  0.99896 

The  situation  is  somewhat  worse  than  this  because  the  velocity  V  is 
determined  at  some  distance  from  the  end  of  the  charge. 

The  effect  just  described  introduces  some  uncertainty  in  cl.assi 
fying  records.  For  example,  when  L/d  is  4  or  less,  the  effect 
described  above  should  result  in  a  detectable  curvature  of  the  re¬ 
cords.  This  curvature  must  be  evaluated  and  compared  with  the  ex¬ 
pected  curvature  for  a  final  classification  of  the  record.  The  ap¬ 
propriate  question  concerns  the  change  in  the  velocity  V  over  the 
observed  distance,  usually  2.5  inches.  At  the  beginning  of  the  re¬ 
cord,  the  detonation  has  reached  the  point  z  =  z^  =  DE  (see  Fig.  4) 
so  that  the  observed  velocity  is  V^  =  F-jD.  At  the  end  of  the  record 
z  =  z2  and  V2  =  F£D,  it  being  assumed  that  the  velocity  D  is  a  con¬ 
stant.  Hence  the  change  in  velocity  is  D(F1  -  F2)  and  the  %  change 
with  respect  to  the  velocity  V1  is 

G  =  100  {F1  -  F2)^F1,  (9) 

due  to  the  geometry  of  the  charge. 

The  following  example  illustrates  the  use  of  the  above  results. 
A  5  inch  diameter  charge  gave  a  velocity  of  4.64  mm/|isec  at  a 
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distance  of  6.7  cm  from  its  free  end.  The  total  length  of  the  charge 
was  10  Inches  or  25  cm.  So  at  6.7  cm  from  the  end,  L/d  is  about  2.4 
and  «  0.978.  At  the  end  of  the  charge  where  L/d  ■  3.3,  Fg  « 

0.9885  and  G  =  -1.1,  which  means  that  the  velocity  decreases  about 
1%  due  to  the  effect  of  geometry.  Analysis  of  the  record  showed  that 
the  velocity  V  decreased  about  256  during  the  observation.  Thus  the 
observed  decrease  is  greater  than  the  expected  decrease  so  that  the 
detonation  could,  in  principle,  be  failing. 

The  correction  described  above  can  be  expedited  somewhat  by  use 
of  the  curves  shown  in  Fig.  5.  The  curves  are  for  25  cm  long  charges 
for  which  the  last  X  cm  are  observed  by  the  streak  camera.  The 
abscissa  in  the  Fig.  5  Is  X  and  the  ordinate  is  -G,  the  %  change  in 
velocity  expected  from  the  geometry  effect. 

Most  of  the  model  propellant  mixtures  being  investigated  have 
been  boostered  by  a  charge  of  high  explosive  such  as  tetryl  or  pento- 
lite.  The  detonation  velocity  of  the  booster  is  greater  than  that  of 
the  acceptor  and  than  that  given  above.  The  derivation  of  the  cor¬ 
rection  is  based  on  the  assumptions  that  the  detonation  front  is 
spherical  not  only  while  it  is  in  the  booster  but  also  after  it  is 
entirely  in  the  acceptor  charge  and  that  the  detonation  velocity  in 
the  acceptor  is  constant.  It  is  recognized  that  these  assumptions 
over-simplify  the  situation.  The  situation  is  diagrammed  in  Fig.  6 
which  shows  the  booster  as  having  the  same  diameter  as  the  test  charge. 
The  length  of  the  booster,  DE  Is  represented  by  h,  the  radius  of  the 
charges  by  r  and  the  detonation  front  by  the  arc  BCD.  Assume  that 
the  center  of  the  circle  1b  at  the  point  0,k),  and  that  its  radius  is 
R.  Substituting  the  coordinates  of  points  C  and  D  Into  the  equation 


of  a  circle  gives 

(H  -  k)2  =  R2 

(10) 

and 

r2  +  (h  -  k)2  =  R2, 

(11) 

where  H  is  the  distance 

AC,  Solving  for  k  gives 

k  =  H2  -  h2  -  r2 

- S(ff— TiF 

(12) 
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so  that/  the  apparent  origin  of  the  spherical  front  is  determined. 

The  quantity  H  is  eliminated  by  noting  that  the  same  intervals 
of  time  are  required  for  the  detonation  to  reach  points  C  and  D.  Let 
.0^  and  be  the  detonation  velocities  in  the  acceptor  and  booster 
charges  respectively.  The  results  are 


V 


h2  +  rg  =  h  +  (H  -  h), 
Di>  Db  Da 


from  which 


where 


H  =  h(l  -  N)  + 


+  r 


(15) 


(14) 


N  =  Da/Db. 

Combining  equations  12  and  14  gives 


k  =  (1-N) 


-  r 


(1  +  N) 


(15) 


2N  Vh2  +  r2  -  h 

Note  that  k  vanishes  when  N  =  1,  i.e.,  when  the  booster  has  the  same 
detonation  velocity  as  the  acceptor.  When  N  <1,  k  is  negative  and 
the  effective  booster  length  is  (h  -  k). 


The  dependency  of  k  on  N  for  1.0,  2.0  and  5.0  inch  diameter 
charges  when  boostered  by  either  1.0  inch  or  2.0  inch  long  booster 
pellets  is  shown  in  Table  2  and  Fig,  7.  The  value  for  k  varies  more 
with  charge  diameter  when  the  booster  length  is  1.0  inch,  than  when 
the  length  is  2.0  inches.  For  N  =  0.5,  the  effective  length,  (h  -  k), 
of  a  2.0  inch  long  booster  is  about  4.2  inches.  Such  a  booster  is 
normally  used  with  an  8.0  inch  long  acceptor  charge,  so  that  the  ef¬ 
fective  charge  length  is  about  12.2  inches.  When  9.0  inch  long 
charges  are  used  with  1.0  inch  long  boosters,  the  effective  lengths 
are  11.1  Inches  for  1.0  inch  diameter  and  11.6  inches  for  5.0  inch 
diameter  charges. 


When  N  <1,  the  effective  charge  length,  Lg,  is  greater  than  the 
actual  charge  length,  see  Table  5.  Thus  the  use  of  the  actual  charge 
length  (donor  plus  acceptor)  in  Eq.  15  results  in  values  of  F  which 
are  too  small.  Use  of  the  effective  charge  length  gives  larger  values 
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of  P  and  results  In  a  more  nearly  correct  value  of  the  axial  velocity. 
The  following  example  may  be  useful  In  demonstrating  these  corrections. 
Assume  the  charge  Is  3»0  inches  in  diameter,  and  that  the  donor  is  1.0 
inch  long.  Thus  for  N  =  1,  the 


TABLE  2 

k  in  inches  for  1.0  and  2.0  inch  long  boosters 


1.0  inch  long  booster  2 

! .0  inch  long  booster 

CHARGE 

DIAMETER. 

inch}  1#0 

mm  EPV 

WmMmm  fcmmm 

1.0 

2.0 

3.0 

N 

k 

k  k 

k 

k 

k 

0.2 

-4.28 

-4.99  -5.93 

•8.15 

-8.57 

-9.20 

0.4 

-1.62 

-1.93  -2.34 

•3.06 

-3.25 

-3.52 

0.6 

-0.73 

-O.89  -1.09 

■1.37 

-1.46 

-1.60 

0.8 

-0.27 

-0.34  -0.43 

■0.51 

-0.55 

-0.61 

1.0 

0.00 

0.00  0.00 

0.00 

0.00 

0.00 

1.2 

0.19 

0.24  0.31 

0.34 

0.38 

0.425 

1.4 

0.33 

0.43  0.56 

0.59 

0.65 

0.74 

effective 

length,  Lg,  is  the  same  as  the  actual 

combined  length  of  the 

charges,  ; 

and  F  =  0.9887.  For  smaller  values  of 

N,  L 

is  increased. 

and  P  is  : 

Increased,  that  is 

,  the  observed  velocity  is 

more  nearly  equal 

to  the  axial  velocity. 

TABLE  3 

Correction  to  Observed  Velocity 

N  k 

L  L  /d 

P 

e  e 

inches 

1.0  0 

10  3.33 

0.9887 

0.8  -0.4 

10.4  3.47 

0.9896 

0.6  -1.1 

11.1  3.70 

0.9909 

0.4  -2.3 

12.3  4.10 

0.9926 

0.2  -5.9 

15.9  5.30 

0.9956 

The  ' 

observed  velocity 

is  about  1; %  too  great 

for  a 

9.0  inch  long. 

5,0  inch  i 

diameter  charge  of 

ammonium  perchlorate  having  a  detonation 

velocity 

of  about  4  mm/nsec 

when  boostered  by  tetryl  having  a 
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detonation  velocity  of  7.2  mm/psec. 

V.  BIFi'EREMTIATING  THE  STREAK  CAMERA  RECORD.  After  a  record  passes 
a  visual  screening  for  lack  of  curvature  it  must  be  differentiated. 

This  usual. y  requires  that  the  record  be  reduced  to  digital  data.  If 
a  graphical  method  is  used  the  data  are  plotted  and  a  curve  is  drawn 
through  the  points.  Next,  normals  or  tangents  are  constructed  at 
intervals  along  the  curve  arid  their  slopes  are  determined.  Thus  the 
derivative  is  determined  at  as  many  points  along  the  curve  as  desired. 

Graphical  differentation  has  been  discarded  because  it  is  a  time 
consuming  process  and  the  results  are  often  unsatisfactory.  There  is 
a  degree  of  subjectivity  in  drawing  a  continuous  curve  through  the 
plotted  points  which  makes  reproducibility  difficult.  Constructing 
the  tangents  (or  normals)  is  also  subjective  and  Is  a  source  of  error. 
However,  graphical  methods  may  be  useful  especially  during  the  early 
part  of  an  investigation  when  the  behavior  of  the  derivative  is  un¬ 
known.  If  the  curve  drawing  is  done  both  carefully  and  honestly, 
abrupt  changes  may  be  detected  in  the  slope  of  the  curve.  These  abrupt 

changes  may  be  of  considerable  interest,  as  in  the  case  of  attenuation 
A 

experiments  . 


Graphical  methods  depend  on  the  accurate  measurement  of  angles. 
One  of  the  difficulties  with  measuring  angles  is  illustrated  by  the 
following.  Assume  the  angle  9  is  measured.  Fig.  8.  The  velocity  is 
proportional  to  k  tan  9,  when  k  is  a  constant  of  proportionality  which 
is  assumed  to  be  well  known.  The  relative  error  in  $  is 


R  =  100  dv  =  100  0 

v  sin#'  cosfl 


200  d0 


(16) 


which  has  a  minimum  at  9  =  45°  where  it  Is  twice  the  error  in  0.  Thus 
the  best  determinations  of  velocity  are  possible  when  9  =  45°.  The 
behavior  of  R/d9  is  shown  in  Fig.  9.  Records  for  which  30°  <  0  <  60° 


are  probably  suitable  for  data  reduction  by  measuring  the  angle.  The 
problem  is  discussed  in  more  detail  by  A.  S.  Bubovik^. 


The  availability  of  electronic  digital  computers  has  increased 
the  use  of  polynomials  and  other  functions  for  fitting  data.  Poly¬ 
nomials  are  easily  differentiated  and  a  table  of  the  values  of  the 
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derivative  is  provided  by  the  computer.  If  the  scale  factors  are 
also  included  in  the  input  to  the  computer,  the  velocity  is  given  as 
a  function  of  one  of  the  variables,  x  or  t,  in  the  output.  Some 
comments  should  be  made  concerning  this  procedure. 

For  one  thing,  the  phenomena  being  studied  may  not  obey  a  poly¬ 
nomial  representation.  In  the  simple  problem  being  discussed  here, 
the  velocity  is  expected  to  approach  a  constant  value  for  reasonable 
values  of  L/d,  A  polynomial  does  not  extrapolate  to  give  a  derivative 
which  approaches  an  asymptate.  Hence  the  choice  of  a  polynomial  may 
be  poor  when  the  physics  of  the  situation  is  considered. 

Polynomials  can  be  made  to  fit  the  data  very  well  simply  by 
using  higher  order  terms.  When  the  degree  is  3  or  greater,  one  or 
more  inflection  points  are  possible.  These  almost  invariably  spoil 
any  extrapolation,  and  they  may  introduce  changes  in  the  derivatives 
which  are  not  expected  from  the  physics  of  the  situation.  Great  care 
should  be  used  in  the  interpretation  of  polynomial  fits  to  data, 
especially  when  one  polynomial  is  fitted  to  the  entire  set  of  data. 

A  piecewise  fit  may  sometimes  give  more  useful  results.  Over  small 
ranges  of  the  data,  even  first  and  second  degree  polynomials  may  be 
useful.  For  the  present  study,  first  and  second  degree  polynomials 
are  useful  in  determining  the  velocity  because  (l)  the  data  are  nearly 
linear  and  (2)  comparison  of  the  fits  with  the  data  give  some  indica¬ 
tions  as  to  both  the  linearity  of  the  data.  Sec.  II,  and  the  scatter 
in  the  data.  Sec.  III. 

Some  curvature  In  the  data  Is  expected  because  of  the  geometrical 
effect  discussed  In  Sec.  IV.  This  should  be  more  evident  In  the 
charges  for  which  L/d  is  small  (4  rather  than  10).  It  would  be  well 
if  this  curvature  could  be  detected  by  the  numerical  analysis  of  the 
data.  In  order  to  do  this,  one  may  use  functions  such  as  the  quadratic 
mentioned  above.  Another  function  which  could  be  used  Is 

x  =  At  +  C(1  -  e'Bt).  (17) 

As  t  becomes  large,  the  first  derivative  approaches  a  constant  value, 

A.  This  is  a  desirable  feature  for  the  present  application  -  at  least 
for  charges  which  have  steady  axial  velocities.  The  function  has  been 
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f 


i 


f* 


C 

used  in  this  application  previously0  but  is  not  used  in  the  present 
data  reduction  operations.  It  could  possibly  be  used  in  place  of  the 
quadratic  for  determining  the  radius  of  curvature  of  a  trace,  see 
Sec.  II. 

It  was  noted  above  that  the  approach  of  the  first  derivative  of 
Eq.  (17)  to  a  constant  value  was  desirable*  Other  functions  have  this 
same  desirable  behavior.  Some  can  be  found  in  the  analysis  of  one 
dimensional  shock  waves,  two  of  which  are  discussed  in  the  following. 
Only  one  of  them  has  been  used  in  the  present  program  of  data  reduction. 

The  decay  of  a  shock  wave  due  to  a  following  rarefaction  wave  is 
discussed  by  Lax^.  The  path  of  the  shock  wave  is  given  by 

At(x-t )+Bt+Cx  =  0  (l8) 

where 


C  =  -k  (uj  +  c-^  -  1) . 


In  the  above  relations,  u^  end  c^  are  the  shock  Induced  particle 
velocity  and  the  shocked  induced  sound  speed,  respectively,  and  k  is 
an  equation  of  state  parameter.  The  interesting  feature  of  this 
representation  is  the  cross  product  of  the  variables.  Means  are 
available  for  fitting  the  constants  to  the  data  in  the  least  square 
sense.  The  function  actually  being  used  is 

x  =  -(A  +  Bt  +  Ct2)  /  (D  +  t).  (19) 

Its  derivatives  are 

dx/dt  =  -C  +  (A  +  CD2  -  DB)  /(D  +  t)2  (20) 

and 

d2x/dt2  =  2 (DB  -  CD2  -  A)  /  (D  +  t)3  (21) 

One  objectionable  feature  of  Eq.  19  is  the  singularity  in  the  function, 
and  in  its  derivatives.  There  is  no  inflection  point,  but  there  is  a 
change  of  sign  in  the  second  derivative  at  the  singularity.  This 
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shows  up  In  some  applications  of  the  function,  i.e.,  the  parameter  D 
is  negative  and  its  absolute  value  is  in  the  range  of  the  values  of 
the  independent  variable  t.  In  the  neighborhood  of  the  singularity 
the  calculated  values  of  the  velocity  and  acceleration  do  not  lie  on 
a  smooth  curve. 

It  is  obvious  that  the  acceleration  (Eq.  21)  vanishes  as  t  becomes 
large.  The  velocity  then  approaches  the  constant,  -C.  Thus  the  func¬ 
tion  behaves  in  the  same  manner  as  the  exponential  function  discussed 
previously.  It  is  this  behavior  which  makes  the  function  of  interest 
in  the  present  context.  Because  it  has  four  parameters,  it  should 
fit  the  data  somewhat  differently  than  Eq.  17. 


Another  derivation  of  the  path  of  a  decaying  shock  is  given  by 

O 

Fowles  .  His  results  are  given  by  the  expression 

C0B  (T-T^^2  +  ACo  (T-T-l)  -  B(X  -  X^  (T  -  T1)1/2  +  (X  -  X^ )  =  0,  (22) 

where 


A  =  -(27  +  3) 

B  =  [oo  -  2(7  +  l)1 


and  C0, 


T-,  are  constants. 
3 


V  y’ 

product  in  common  with  the  Lax  expression, 
plied  in  data  reduction. 


'  T1 

This  representation  has  a  cross 
It  has  not  yet  been  ap- 


The  work  Involved  in  using  such  functions  as  the  hyperbolic  and 

the  exponential  described  above  is  greatly  expedited  by  the  use  of  a 

subroutine  "A  General  Weighted  Least  -  Square-Fit  Subroutine".^  A 

useful  feature  of  the  subroutine  is  the  calculation  of  the  quadratic 

mean  error  (QME)  for  the  overall  fit  and  for  each  of  the  parameters 

Involved.  The  subroutine  can  be  used  in  different  ways  to  determine 

the  value  of  the  parameters  of  equations  such  as  Eq.  18.  If  the 

equation  is  fitted  as  it  is  written,  the  values  of  the  QME  for  some 

of  the  parameters  may  be  enormous.  This  is  caused  by  the  singularity 

which  gives  large  numbers  for  the  matrix  which  must  be  inverted  in  the 

subroutine.  The  trouble  is  circumvented  by  regarding  both  x  and  t  as 

independent  variables,  and  the  cross-product  term,  xt,  as  the  dependent 

variable. *  This  approach  produces  more  reasonable  values. for  both  the 
*This  method  was  suggested  by  J.  W.  Forbes. 
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parameters  and  for  their  QME's.  The  QME  of  the  fit  is  still  calcu¬ 
lated  on  the  basis  that  t  is  the  independent  variable  and  that  x  is 
the  dependent  variable. 

Another  method  for  handling  the  data  uses  a  spline  function. 

This  consists  of  a  set  of  cubic  functions  which  are  fitted  to  the 
data  in  a  least  squares  sense.  Each  cubic  represents  a  subset  of  the 
data,  and  the  first  and  second  derivatives  are  required  to  be  continu¬ 
ous  at  the  joints  between  adjacent  cubics.  The  set  of  data  may  be 
divided  into  subsets  of  3  or  more  points.  When  only  a  few  points  are 
used  in  a  subset,  the  spline  function  employs  many  parameters  so  that 
local  variations  in  the  data  may  be  followed.  This  ability  to  fit  the 
data  closely  is  a  characteristic  which  the  spline  shares  vrith  higher 
degree  polynomials.  The  spline  function  smooths  the  data  in  about  the 
same  sense  as  the  usual  smoothing  formulas  frequently  used  for  data 
whi  h  are  equally  spaced  in  one  variable.  Stated  another  way.  the 
spline  function  permits  the  smoothing  of  data  in  those  situations 
where  it  is  thought  inadvisable  or  impossible  to  advance  one  variable 
in  equal  steps.  The  spline  function  parameters  are  determined  by  a 
computer  program  which  also  evaluates  the  derivative.  Thus  the  pro¬ 
gram  accomplishes  the  same  result  as  the  more  usual  smootning  and 
differentiation  operations  on  equally  spaced  data. 

VI.  RESULTS  OF  NUMERICAL  DATA  REDUCTION.  Some  of  the  operations 
which  may  be  performed  on  the  data  were  described  in  the  previous 
sections.  These  operations  are  performed  by  a  digital  computer  and 
the  results  are  displayed  in  tabular  form  and  in  graphical  form.  The 
table  and  the  plots  will  be  presented  and  described  in  this  section. 
The  meaning  of  the  results  will  be  discussed  in  the  following  section. 

The  first  line  of  the  computer  output.  Table  4,  gives  the  shot 
number  to  identify  the  data.  The  second  line  gives  conversion  factors 
so  that  measurements  from  the  Telereader  can  be  converted  to  physical 
units.  The  raw  data  are  given  in  columns  3  and  4  of  the  columnar  part 
of  the  table.  The  converted  data  are  given  in  column  1  and  2,  the 
units  being  psec  for  time  and  mm  for  distance. 
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The  crude  values  for  the  lnstanteous  velocities  are  given  in  the 
column  headed  DXDT.  Values  in  this  column  should  vary  very  little  if 
the  streak  camera  record  is  straight  and  easy  to  read.  Hence  any 
value  of  DXDT  differing  significantly  from  the  average  signals  a  bad 
data  point. 

The  lnstanteous  velocity  is  given  in  columns  5  and  7  also.  The 
column  labeled  SPLINE  gives  the  velocity  as  evaluated  from  the  spline 
function,  (SF).  The  column  labeled  HYPERBOLIC  gives  the  velocity  as 
calculated  from  the  hyperbolic  function  (HF)  Eq.  19.  The  acceleration 
as  given  by  the  HF  is  given  in  column  8.  Residuals  for  the  three 
functions  are  given  in  the  last  3  columns  of  the  table. 

Coefficients  for  the  linear,  quadratic  and  HF  are  given  near  the 
top  of  the  table  (XT  FIT  means  the  hyperbolic  function).  Immediately 
below  the  value  of  each  coefficient  is  its  percentage  error  based  on 
the  quadratic  mean  error  (QME)  of  that  parameter  which  Is  essentially 
the  standard  deviation.10  Where  the  error  in  a  parameter  is  large — 
100$  or  more — the  term  Involving  that  parameter  should  probably  be 
discarded  and  a  new  fit  should  be  obtained. 

Below  the  columnar  part  of  Table  4  are  results  which  help  to 
determine  if  the  record  represents  a  steady  detonation.  The  parameters 
of  the  quadratic  are  used  to  give  the  acceleration  and  the  time  and 
distance  required  for  a  10$  change  in  the  velocity.  On  the  next  line 
is  the  percent  change  in  the  velocity  over  the  time  of  observation, 
i.e.,  the  Interval  of  time  between  the  first  and  last  entries  in 
column  1, 

The  root  mean  square  of  the  error  in  the  spline  function,  (SF), 
is  given  in  the  third  line  from  the  bottom  of  the  table.  The  deriva¬ 
tive  of  the  SF  can  be  obtained  by  the  use  of  a  subroutine  for  any  value 
of  the  time.  For  example,  the  subroutine  was  used  to  obtain  the  de¬ 
rivative  (the  velocity)  at  each  value  of  the  time  in  column  1  of  the 
table.  These  values  are  recorded  in  column  5.  Another  set  of  values 
are  obtained  at  45  equally  spaced  values  of  the  time  within  the  range 
of  the  entries  of  column  1.  These  results  are  divided  into  three 
groups  of  15  each  and  the  average  Is  calculated  for  each  group  as  well 
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as  for  the  set  of  45  values.  These  averages  are  reported  on  the 
second  line  from  the  bottom  of  the  table.  The  difference,  in  percent, 
of  the  velocity  over  each  one-third  with  respect  to  the  overall  aver¬ 
age  is  given  on  the  last  line. 

As  an  aid  in  examining  the  results,  plots  are  prepared  showing 
DXDT  and  the  velocity  as  computed  from  the  spline  function,  see 
Pig.  10.  In  the  plots,  values  of  DXDT  are  represented  by  points, 
while  the  results  of  the  SP  are  given  as  connected  points. 

VII.  DISCUSSION  OP  RESULTS.  The  results  presented  in  Table  4  and 
Pig.  10  for  Shot  487  will  be  discussed  first.  Later,  the  results  for 
other  shots  will  be  presented  and  discussed. 

Shot  487  gave  a  good  record  as  is  evident  from  examining  the 
values  of  DXDT  in  Table  4  and  in  Pig.  10.  Note  that  in  the  table  the 
first  value  of  DXDT  is  arbitrarily  set  to  zero.  This  results  in  a 
problem  with  the  plotting  procedures  which  is  circumvented  by  roaking 
DXDT  (l)  =  DXDT  (2)  for  the  plot.  The  values  of  DXDT  dc  vary,  having 
a  minimum  value  of  2.65  and  a  maximum  of  5-02.  The  velocity  reported 
for  this  shot  was  2.86  mm/psec,  the  value  of  B  for  the  linear  function. 
Hence  the  maximum  variation  of  DXDT  with  respect  to  the  reported 
velocity  is  +5$  and  -rJ%,  Most  of  the  values  lie  within  see 

Fig.  10.  These  variations  are  tolerable  and  it  is  probably  not 
necessary  to  discard  the  few  points  which  deviate  more  than  5 %• 

For  the  results  given  here,  the  data  were  divided  Into  subsets 
of  four  points  each  in  obtaining  the  fit  by  the  SF»  Used  in  this  way, 
the  SP  smooths  the  data  considerably  so  that  the  velocities  from  the 
SF  vary  less  than  do  the  values  of  DXDT.  The  SP  preserves  the  general 
trend  of  the  velocity,  however.  If  greater  smoothing  had  been  de¬ 
sired,  the  number  of  data  points  in  each  subset  would  have  been  in¬ 
creased.  Note  that  the  values  from  the  spline  function  (SF'j  are  low 
at  the  beginning  of  the  record,  as  are  the  values  of  DXDT.  This  is 
frequently  observed  in  the  data  and  the  cause  is  not  known.  The  im¬ 
plication  is  that  the  detonation  is  accelerating.  This  is  not  ex¬ 
pected,  especially  in  view  of  the  results  presented  in  Sec.  IV,  which 
predicts  a  decreasing  velocity.  Also  note  that  the  velocity  drops 
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near  the  end  of  the  record  so  that  the  average  over  the  last  1/3  is 
lower  than  the  average  over  the  1st  and  2nd  thirds.  This  latter  por¬ 
tion  of  the  data  indicates  that  the  detonation  is  decelerating. 

The  increase  of  velocity  at  the  beginning  of  the  record,  and  its 
decrease  over  the  latter  portion  of  the  record,  as  shown  by  the  SP, 
make  the  data  difficult  to  fit.  Note  that  the  residuals  for  the  linear, 
quadratic  and  HP  fits  have  a  cluster  of  negative  signs  near  the  end  of 
the  record  and  that  the  last  three  residuals  are  positive  for  all 
three  functions.  The  maximum  residual  is  about  0.4#  of  the  value  of 
the  distance.  These  small  values  may  be  meaningless  so  that  one  can, 
for  example,  regard  the  linear  fit  to  be  adequate.  Hence  the  velocity 
of  the  detonation  is  2.86  ±  0.1#  (see  the  value  of  B  for  the  linear 
fit).  However,  the  last  5  values  of  the  velocity  from  the  SP  are  less 
than  2.86  so  it  is  not  clear  if  the  charge  detonated  at  constant 
velocity. 

The  quadratic  gives  a  velocity  of  2.86  mm/usec  at  11  nsec  (the 
middle  of  the  record)  in  agreement  with  the  value  from  the  linear  fit. 
Because  the  value  of  C  is  negative,  the  quadratic  also  shows  that  the 
detonation  velocity  is  decreasing,  in  agreement  with  the  SF.  But  ob¬ 
serve  the  errors  for  B  and  C  for  the  quadratic,  0.5#  and  79#  respec¬ 
tively.  The  latter  is  so  large  that  little  confidence  can  be  placed 
in  the  prediction  of  the  acceleration  by  the  quadratic.  Hence  these 
computations  have  been  of  little  help  in  determining  if  the  detonation 
was  steady. 

The  data  and  results  for  Shot  502  are  given  in  Table  5  and  Fig.  11. 
The  charge  was  5.08  cm  in  diameter  and  about  25  cm  long  so  that  the 
observed  velocity  should  be  decreasing  slightly  (see  Sec.  IV).  This 
shot  is  of  interest  because  the  quadratic  indicates  the  detonation  is 
accelerating  2.3#  over  the  time  of  observation.  The  error  in  the  co¬ 
efficient  C  for  the  quadratic  is  25#  here  compared  to  79^  for  Shot  487 
which  was  discussed  previously.  Hence  the  quadratic  fits  the  data  of 
Shot  502  better  than  that  of  Shot  487,  and  the  calculated  acceleration 
is  more  reliable. 
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A  better  fit  to  the  data  of  Shot  802  could  possibly  be  obtained 
by  discarding  the  data  on  the  lines  where  the  values  of  DXDT  are  about 
3.47  and  3.52. 

The  velocities  for  t  =  10  jxsec  are:  (l)  linear,  3.304,  (2) 
quadratic,  3.305*  (3)  HP,  3.290,  (4)  SF,  3.365.  The  agreement  among 
the  values  is  good,  and  is  improved  if  the  average  over  the  2nd  third 
of  the  record  is  used  for  the  SP,  i.e.,  3.307.  Note  that  the  results 
obtained  by  averaging  the  results  from  tne  SP  show  acceleration.  The 
"ripple"  which  is  in  evidence  in  Pig.  10  is  an  interesting  feature  of 
this  record. 

Data  and  results  for  Shot  504  are  given  in  Table  6  and  Fig.  12. 
The  charge  was  2.54  cm  in  diameter  and  25  cm  long  so  that  the  observed 
velocity  should  be  close  to  the  axial  velocity.  The  record  was  dif¬ 
ficult  to  read — see  the  variations  in  the  values  of  DXDT.  The  veloci¬ 
ties  are:  (1)  linear,  3.872  ±  0.2g,  (2)  quadratic,  5.873  A  0.7?', 

(3)  SF,  6.027,  (4)  HP,  5.843  where  the  last  3  values  are  taken  at  t  -- 
5.25  usee.  Agreement  is  remarkable,  except  for  the  value  from  the  SP. 
Averaging  over  the  middle  third  of  the  record  gives  5.999,  which  is 
still  larger  than  any  of  the  other  results. 

The  error  in  the  acceleration,  i.e.,  i  C  for  the  quadratic  fit, 
is  so  large  that  the  quadratic  is  useless  for  the  prediction  of 
steadiness.  Eliminating  a  few  points  at  the  start  and  at  the  end  of 
the  record  might  permit  a  better  determination  of  the  acceleration. 
Removal  of  these  data  would  eliminate  the  large  values  of  DXDT  at  each 
end  of  the  plot,  Pig,  12. 

Each  of  the  three  records  discussed  contains  a  certain  amount  of 
"ripple".  This  shows  up  both  in  the  plot  of  DXDT  and  in  the  plot  from 
the  SP.  This  ripple  could  be  the  result  of  faulty  bearings  in  the 
streak  camera.  Or  the  film  reader  could  be  defective.  It  is  possible 
of  course,  that  the  detonations  are  really  not  steady.  The  idea  that 
detonations  are  unsteady  is  somewhat  abhorrent  so  that  the  other 
possibilities  are  being  investigated  first. 

In  order  to  check  the  film  reader,  objects  with  straight  edges 
were  placed  In  the  projector.  In  all  cases,  some  ripple  was  observed 
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on  the  plot  of  the  velocity  -  which  was  simply  the  slope  of  the  line 
as  seen  on  the  screen  of  the  machine.  Plots  from  two  of  these  experi¬ 
ments  are  shown  in  Pig.  13.  The  top  plot  is  the  result  obtained  with 
a  stretched  thread  in  the  projector.  There  is  a  ripple  in  both  the 
results  from  the  SP  and  DXDT.  Similarly,  a  ripple  is  observed  for  a 
ruler,  see  the  lower  plot  in  Pig.  13.  These  results  show  that  either 
the  film  reader  is  defective,  or  that  it  is  not  being  used  properly. 

The  trouble  is  thought  to  be  caused  by  the  illuminating  system  in  the 
machine  which  consists  of  several  small  fluorescent  lamps.  There  is 
no  diffuser  or  condensing  lenses  between  the  lamp  and  the  film.  The 
screen  illumination  is  nonuniform  to  such  a  degree  that  the  operator 
may  find  it  impossible  to  set  the  cross  wires  on  the  image  properly. 

The  machine  will  be  replaced  in  the  near  future  and  improved  results 
may  then  be  obtained. 

In  some  cases  the  errors  in  the  coefficients  of  the  hyperbolic 
function  are  large.  However,  the  residuals  may  be  as  small  as  those 
for  the  linear  and  quadratic  functions.  The  meaning  of  the  large 
uncertainties  is  that  the  terms  involving  A  and  D  could  be  dropped. 

The  result  would  then  be  x  +  B  +  Ct  =  0.  That  is,  a  linear  function 
is  more  appropriate  for  the  data.  It  is  of  interest  to  examine  the 
behavior  of  the  hyperbolic  function  when  non-linear  data  are  used. 
Results  were  reported'5  for  5.08  cm  diameter  tetryl  charges  in  work 
concerned  with  the  correction  for  the  finite  length  of  charges.  Data 
for  three  of  the  experiments  have  been  treated  in  the  same  way  as  the 
propellant  data  (see  Tables  7,  8,  and  9).  Because  the  original  data, 
where  the  units  were  counts,  were  not  available,  the  previously  con¬ 
verted  data  were  used  as  input.  Hence  the  3rd  and  4th  column  are  the 
same  as  the  1st  and  2nd  column  respectively  and  distance  is  in  cm 
rather  than  mm.  There  is  at  least  one  bad  pair  of  values  of  x  and  t 
in  Table  7  (see  the  value  1.086  in  the  DXDT  column).  This  bad  data 
set  is  not  smoothed  very  well  by  the  SF.  Curvature  is  so  great  that 
the  linear  function  is  useless.  The  HF  fits  the  data  fairly  well  (see 
the  residuals  for  the  three  tetryl  charges,  Tables  7,  8,  and  9  and 
the  values  for  the  QME  in  Table  10.)  The  coefficients  and  their  errors 
are  collected  in  Table  10,  as  well  as  results  for  the  three  propellant 
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shots  discussed  previously.  Smaller  errors  are  observed  in  A  and  D 
for  the  tetryl  charges  than  for  the  propellant  charges.  That  is, 
these  coefficients  are  more  meaningful  for  the  tetryl  data,  due  to  the 
greater  curvature  of  the  (x,  t)  curve  being  fitted. 

The  average  velocity  from  6  tetryl  charges  was  reported  to  be 
7.17  mm/psec.-^  The  purpose  of  the  referenced  work  was  to  demonstrate 
that  the  steady  detonation  velocity  could  be  determined  by  observing 
the  detonation  of  relatively  short  charges.  Interpretation  of  the 
experimental  data  required  the  use  of  the  function 

x2  =  D2  t2  +  A  D  t  (23) 

in  which  D  and  A  are  parameters.  When  extrapolated  to  large  values  of 
x  and  t,  this  function  gives  the  steady  velocity  D.  The  same  deter¬ 
mination  can  be  made  with  the  HF  being  discussed  here.  The  velocity 
in  question  is  simply  the  value  of  -C  for  the  HF.  For  tetryl  shots  1, 
2  and  3  the  velocities  are  7*^9,  8.52  and  f.72  mm/psec  respectively. 
These  values  are  all  greater  than  those  reported  previously,  so  that 
Eq.  23  gives  better  results  than  the  HF.  This  is  not  particularly 
disappointing  because  Eq.  23  is  based  on  physically  reasonable  assump¬ 
tions  for  the  particular  situation,  i.e.,  charges  of  finite  dimension. 
The  HF  is  also  based  on  physically  reasonable  assumptions,  but  the 
situation  being  described  involves  strictly  one -dimensional  flow. 

VIII.  CONCLUSIONS. 

For  a  good  streak  camera  record,  the  velocity  can  be  obtained  in 
any  of  several  ways.  That  is,  the  results  from  the  linear,  quadratic, 
hyperbolic  and  spline  functions  all  agree  reasonably  well.  For  the 
spline  function,  it  is  necessary  to  average  over  all,  or  a  portion,  of 
the  results.  It  Is  possible  that  averaging  could  be  avoided  by  stif¬ 
fening  the  spline  function,  i.e.,  the  data  can  be  treated  in  larger 
subsets  so  that  the  spline  function  would  take  on  the  characteristics 
of  a  cubic  function. 

One  of  the  main  needs  az  present  is  a  means  of  determining  if  the 
velocity  is  steady.  The  best  means  found  so  far  are  the  crude  in¬ 
stantaneous  velocities,  the  column  DXDT  in  the  computer  output,  and 
the  velocities  from  the  spline  function.  However,  the  ripple  ir  the 
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velocity  makes  the  problem  difficult.  Replacement  of  the  Telereader 
may  make  this  task  easier.  If  the  present  Telereader  is  introducing 
the  ripple  into  the  data,  the  accelerations  computed  from  both  the 
quadratic  and  the  hyperbolic  function  may  become  useful. 

One  way  of  using  the  results  of  the  spline  function  for  testing 
the  steadiness  of  the  velocity  is  to  examine  the  relation  between  the 
overall  average  of  the  velocities,  and  the  averages  over  each  one 
third  of  the  results.  This  comparison  suffers  from  the  ripple  also. 
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